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МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ 
СТРАТИФИЦИРОВАННОГО ТЕЧЕНИЯ КРОВИ 

В КРОВЕНОСНОМ СОСУДЕ В ВИДЕ МОДЕЛИ 
РАЗДЕЛЬНОГО ТЕЧЕНИЯ ДВУХ НЕСМЕШИВАЮЩИХСЯ 
ЖИДКОСТЕЙ В СУЖАЮЩЕЙСЯ КОНИЧЕСКОЙ ТРУБКЕ 


В работе на основе модели раздельного течения двух несмешивающихся жидкостей 
в сужающейся конической трубке решена задача о математическом моделировании 
стратифицированного течения крови в кровеносном сосуде. Дана оценка влияния 
вязкостных отношений слоев на профили распределения скоростей и на относи- 
тельный расход. 
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стенчатый эффект. 


Как известно [1-3], кровь, движущаяся по кровеносным сосудам, 
может рассматриваться как суспензия (жидкость со взвешанными в ней ча- 
стицами). Течение такой суспензии обладает такими свойствами, что в уз- 
кой зоне около стенок сосуда взвешанные частицы отсутствуют. Это явле- 
ние носит название пристенчатого эффекта. При этом концентрация жид- 
кости большой вязкости практически равна нулю у стенок сосуда и макси- 
мальна в окрестности ее оси. Таким образом, учет пристенчатого эффекта 
требует при изучении движения крови в кровеносных сосудах учесть слои- 
стый характер ее движения, опираясь на уравнения Навье-Стокса. 

В данной работе решается задача о стратифицированном (слои- 
стом) течении двух несмешивающихся жидкостей в сужающейся кониче- 
ской трубке (рис.1). 





Рис.1. Схематическое изображение слоистой жидкости в конической трубке 


В качестве ИСХОДНЫХ уравнений берется линейная система уравне- 
ний Навье-Стокса для случая «тонкого слоя» для ВЯЗКОЙ несжимаемой жид- 
КОСТИ: 
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где ®.,0, - компоненты вектора скорости; и; - динамический коэффици- 


ент вязкости; р; — давление в слоях; г 2 - цилиндрическая система 
коэффициентов. 
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Система уравнений (1) решается при следующих граничных услови- 
ях: 
— условие прилипания жидкости к стенкам трубки; 


— ограниченность функции 9, ид, на оси трубки; 


— на границе раздела слоев равенство нормальных и касательных 
напряжений и условие (стратифицирование) раздельного течения жидко- 
сти. 

Входное и выходное давление считается заданными. 


Введем функцию тока И (7.2) по формулам: 


ИЕ (2) 
92 ду 


С учетом (2) уравнения (1) можно записать в виде: 
9’? 19у, 19 19, 1 р 
Е ва м (3) 


— = 
о! и, 92 

Точное автомодельное решение задачи. Решение системы уравнений (3) 

будем искать в виде 


1 ар. тв г 
и, И, ©), Е 6 = -а2‘ (4) 
и, 42 е е 


1 


Подставив (4) в (3), получим: 
Та, 1 тай. 
к 14, 14а 14а, = с (= 1,2). (5) 

7% 6% о 4% 
Граничные условия при этом записываем в виде: 
4, ау, _ ау 
= Опри & = 4; = 2 при ё = 4 
Фу, ау, ау 
= приё = А; — "= 0 приё = 0. (6) 
И 4 2 Д 2 4 2 р 6 1 ё р ) 


-( -( 
Здесь у= А.е “, г= Ае “” - соответственно уравнения 


























контура стенки конической трубки и границы раздела слоев жидкости 
(4> А], а = с0п54. 


Условия раздельного (стратифицированного) движения на границе 
раздела слоев автоматически выполняются, поскольку 


| -(2 
—. = -Ае““а приё = А, 
9, 


поэтому на границе раздела вектор скорости направлен по касательной к 
контуру раздела слоев жидкостей. 

Решение задачи (5)-(6) находим непосредственным интегрировани- 
ем, в результате получаем: 
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3 
Е + СЕ ШЕЕ. (1) 





в = ра Норыё чей, В, = 


а. граничные условия (5) ДЛЯ ое. ПОСТОЯННЫХ ИН- 
тегрирования, придем к следующей системе алгебраических уравнений: 








3 
- + САША, + СА, = 0, 





3 3 
сА с›А 
с = 0, А+ с,А = — 1 +сАША+сА, 
4 
3 3 
Ао + С. = Ё2 -с, А+ с ША с. (8) 
1 
Решив систему (8), получим: 
Ш. Р ых Рех . 
И -- Я @2\? 
| и ›(е'-е°) 
2 2 2 
А 4 _ с.Аг. 


2 


= 0: ВЫ ов. о 





2 я 
А 1. А, Ш | 
4 4 ° с. А? 
- = с,шА,. 9) 
И > Ш + й > - Ш 4 
И 1 А, И 1 А, 
Расход в любом сечении трубки определяется выражением: 


О= О, + О,, (10) 


С 


с; = 





где 
А 4 2 

С уз с. , с, 

ЕО ОС о 

ото Чрн -Чь 
А) 2 2 
А А 
- 2+ сё ШС а = А - Ас — 6—2 
о броне 4-е 44 


2 2 
654 2 2 А, А 

А’ - А+ с. —Ш—; А> А. 
$ 1 ; 57 А 2 1 


С учетом формул (2) и (7) для компоненты скорости 1, получим 





следующие выражения: 
2 
о 1 56” 
й ВИЙ -242 
е 4 


ре А; 
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1 а 
2 
ие Е с —1 с. шё ТС А Ё А.. (11) 
е 4 
Из формул (11) с учетом (9) будем иметь: 
ри А 4 Ав, А 
0 = 4 -- р Чо М ОЕ д; 
2 2 2 2 > 
ВЫ 7 А_4 Ан, А 
4 4 4 и, А, 
2 2 А? 
(2) Ея а + ыы а ше 
0, 4 4 4 А, 
(1) - ы А? А? А? А 2 А, ё А,, (12) 
ово С4 п _ 22 _ Аа 
4 4 4, 4 
ри, АИ, | 
4 
где Л = р й 





о те. 


| 1 А, | 1 А 
Из зависимостей, приведенных на рис.2, следует, что чем меньше 
вязкость слоя жидкости, прилегающая к стенкам сосуда, тем более макси- 
мальное значение относительной осевой составляющей скорости этого 
слоя. 





С учетом формул (11) для относительного значения расхода 


[8 ы О, + О,) / О получим следующее выражение: 

















А АА А 4 А, 
Е 
©, 4-А | А 4, А 4 : № [42-4] 

16 Г: О Е: 


Значения относительного расхода О при различных значениях отношения 
вязкости слоев и от отношения протяженности слоев 
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Из зависимостей, приведенных в таблице, следует, что чем меньше 
| — ‚ тем меньше относительное значение расхода 0. При этом с уве- 
1 


личением толщины пристенчатого слоя относительное значение расхода 
[2 резко возрастает. Как и следовало ожидать, при Ш. = Ц; в предель- 


ном случае 4 > 0, О 0, О- 1. 


кз! а 


—__ 
м 
вы ы а к. 
в. ба. > 7? 
0.04 ьч д. = = 
м рии 
и. == © 


0.02 
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Рис. 2. Эпюра распределения относительной осевой 
составляющей скорости (А = 0,03см, 4-А= 0,04 см) 
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